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We present a measurement of the branching fraction for the rare decays B — ► pev and extract a 
value for the magnitude of V u t, one of the smallest elements of the Cabibbo-Kobayashi-Maskawa 
quark-mixing matrix. The results are given for five different calculations of form factors used to 
parametrize the hadronic current in semileptonic decays. Using a sample of 55 million BB meson 
pairs recorded with the BABAR detector at the PEP-II e + e~ storage ring, we obtain B(B° — ► 
p- e+u) = (3.29 ± 0.42 ± 0.47 ± 0.60) x 10~ 4 and \V ub \ = (3.64 ± 0.22 ± 0.25 tojg) x 10~ 3 , 
where the uncertainties are statistical, systematic, and theoretical, respectively. 

PACS numbers: 13.20.He, 12.15.Hh, 14.40.Nd 



Exclusive b — > ulv decays can be used to determine \V u b\, one of the smallest and least well-determined ele- 
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merits of the Cabibbo-Kobayashi-Maskawa quark-mixing 
matrix. The modes B — > pev have a comparatively 
large branching fraction, and a high fraction of events 
is found at large electron momenta. We determine both 
the branching fraction B(B — > pev) and \V u b\ using form 
factors, which describe the hadronic current in the decay, 
to extrapolate the decay rates to the full range of lepton 
energies and to normalize B to | V u b \ . Five different form- 
factor calculations are used, as given in Table ID 

The data in this analysis were collected with the 
BABAR detector at the PEP-II Q asymmetric-energy 
e + e~ storage ring. The integrated luminosity of the 
sample recorded on the T(4S) resonance in years 2000 
and 2001 ( "on-resonance" ) is 50.5 fb _1 , corresponding to 
55.2 million BB meson pairs. An additional 8.7fb _1 
of data were taken 40MeV below the resonance ("off- 
resonance" ) . BABAR is a detector optimized for the asym- 
metric beam configuration at PEP-II. Charged-particle 
momenta are measured in a tracking system consisting 
of a 5-layer, double-sided silicon vertex tracker (SVT) 
and a 40-layer drift chamber (DCH) filled with a mixture 
of helium and isobutane, both operating in a 1.5-T su- 
perconducting solenoid. The electromagnetic calorimeter 
(EMC) consists of 6580 CsI(Tl) crystals arranged in bar- 
rel and forward endcap subdetectors. Particle identifica- 
tion is performed by combining information from ioniza- 
tion measurements in the SVT and DCH, energy deposits 
in the EMC, and the angle and number of Cherenkov 
photons measured by the DIRC (detector of internally 
reflected Cherenkov light). 

We select decays in the modes B + — > p°e + v, B° —> 
p~e + v, B + — > u>e + v, B + -^n°e + v, and B° — > ir~e + v , 
with p° — > 7r + 7r~, p~ — > ir°ir~ , and ui — > 7r 7r + 7r~. The 
inclusion of charge conjugate decays is implied through- 
out. The analysis is optimized for B — ► pev decays, 
similar to that in Ref. [8j. Each signal event is some- 
times reconstructed in one of the four other modes; the 
7r and lu modes are included in order to estimate this 
crossfeed into the p modes. Throughout this paper, 
all variables are expressed in the T(4S) center-of-mass 
frame, except if stated otherwise. Two electron-energy 
regions are considered: 2.0 < E e < 2.3 GeV (low-E e ) 
and 2.3 < E e < 2.7 GeV (high-E e ). A large background 
to b — > uev decays comes from the more copious b — > cev 
decays. This background is kinematically suppressed in 
the high-i? e region and dominates in the low--E e region. 
The low--E e region provides the background normaliza- 
tion in the high-i? e region. The largest background in 
the high-i? e region is continuum e + e~ — » qq events. The 
off-resonance data are used to estimate its size. 

Hadronic events are selected based on track and pho- 
ton multiplicity and event topology. We use tracks orig- 
inating from the interaction point with at least 12 hits 
in the DCH and a transverse momentum greater than 
0.1 GeV/c. Signals in the EMC with E lab > 30 MeV 
that are not associated with any track are considered as 



TABLE I: Form-factor calculations used in the determination 
of B(B — » pev) and \V u b\, and predicted normalizations Ith 
(as defined later in Eq.^J. 



Form factors 


r th (ps- 1 ) 


Error (%) Reference 


ISGW2 


14.2 


±50 


1 


Beyer/Melikhov 


16.0 


±15 


[2] 


UKQCD 


16.5 


±21,-14 


[3] 


LCSR 


16.9 


±32 


f4] 


Ligeti/Wise 


19.4 


±29 


[5J 



photons if the lateral moment of the shower energy dis- 
tribution is smaller than 0.8. We select events with 
at least five tracks, or with at least four tracks and at 
least five photons. We require the ratio H2/H0 of Fox- 
Wolfram moments to be less than 0.4. This require- 
ment keeps 85% of the pev signal; it rejects 55% of the 
non-BB events. 

Electrons are identified with a likelihood estimator us- 
ing information from the DCH, EMC, and DIRC subde- 
tectors [ll|. The selection efficiency is around 90%, with 
a pion misidentification rate of less than 0.1%. We reject 
electrons from J/ip decays and from photon conversions. 

Charged pion candidates are tracks not identified as 
kaons with high confidence based on DIRC and dE/dx 
measurements. A 7r° is reconstructed from photon pairs 
with an invariant mass 120 < M 77 < 145 MeV/c 2 . 

To reconstruct p° mesons, we combine two oppositely- 
charged pions, and for p^ a pion track and a ir°. To sup- 
press combinatorial background we require that the pion 
with the higher momentum satisfies p^ > 400 MeV/c 
and the other pion p„ > 200 MeV/c. For the u>, we com- 
bine two oppositely-charged pions with a 7r° . To suppress 
combinatorial background we require p^ > 100 MeV/c 
for each pion. In the mode B — > nev we require 
p v > 200 MeV/c. 

The missing momentum in the event is given by 

Pmiss = - ^2 Pi - ^2 Pi , (1) 

tracks photons 

where the sums are over all accepted tracks and photons. 
We require |cos6> m i ss | < 0.9, where 6> m i ss is the angle 
between p m i ss and the beam axis. This rejects events 
with missing high-momentum particles close to the beam 
axis. We also compare the direction of p m i S s with that 
of the neutrino inferred from p v = ps — Py, where Y is 
the p ± e, lo ± e, or n ± e system. The latter is known 
to within an azimuthal ambiguity about the B direction 
since only the magnitude of pb is known. We use the 
smallest possible angle A^min between the two directions 
and require cos A6> m i n > 0.8. Using the constraints Eb — 
Ebcam and pi = {ps — Py) 2 = 0, the angle between the 
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B meson and the Y system is 

2E B E Y -(M*+M*)<t 

cos & BY = ,^ ,73 ; 2 . (2) 

Signal events fulfill |cos0gy| < 1; allowing for detector 
resolution we require |cos6*sy| < 1.1. After all other 
selection criteria, this requirement rejects more than 60% 
of the b — > cev and approximately 68% of the remaining 
continuum backgrounds, it retains 98% of the signal. 

To further reduce the continuum background, we use 
a neural net with 14 event-shape variables: the sum of 
track and photon energies in nine cones centered on the 
lepton-momentum; the angle ^thrust between the thrust 
axis of the Y system and the thrust axis of the rest of 
the event (the thrust axis is defined to be the direction 
that maximizes the sum of the longitudinal momenta of 
all particles); the angle #thrust,y between the thrust of 
the Y system and the beam axis; the angle #i ep t,rest be- 
tween the direction of the lepton and the direction of the 
total momentum of all tracks except the Y system; the 
momentum of the track with the smallest opening angle 
with respect to the electron; YliPi ' ?W Wi\> where n e 
is the direction of the electron and pi are the momenta 
of all tracks except the electron. After all other selection 
criteria, the neural net condition removes more than 90% 
of the continuum events in the high-i5 e region, while re- 
taining approximately 60% of the signal events in each 
signal mode. 

After all selections, there remain on average 3.4 candi- 
dates per event. We choose the one with a total momen- 
tum \py + Pmiss | closest to the i?-meson momentum \p B \ ■ 
The probability of making the right choice for the signal 
modes is approximately 85%. 

The total efficiency in the high-_E e region is 12.0% 
(9.5%) for the mode B+ -> p°e + v (B° -> p~e + v) in 
the ISGW2 model; it is 4.2% (3.3%), when relating the 
accepted events in the high-_E e region to events with all 
electron energies. 

We perform a binned maximum-likelihood fit to 
the two-dimensional distribution (M^y, AE), where 
M ww ( w ) is the invariant mass of the p (ui) meson and AE 
is the difference between the reconstructed and the ex- 
pected B-meson energy, AE = Shadron + E e + |p miss |c - 
Ebeam- The fit is performed simultaneously for the five 
signal modes in the two E e ranges. For the B — > pev 
modes, the data are divided into 10 x 10 bins over the 
(M w , AE) region 0.25 < < 2.00 GeV/c 2 and 

\AE\ < 2 GeV. For the cu channel, we use five bins in 
the range 702 < M TO < 862 MeV/c 2 and ten bins in 
\AE\ < 2 GeV. For the modes B -> uev, only AE is 
used as a fit variable, also with ten bins. 

In the fit, the likelihood is calculated as product of 
probability distributions for each of the five signal modes, 
for other b — > uev decays, for b — > cev decays, for contin- 
uum events, and for a small contribution due to misiden- 
tified electrons. Shapes and normalizations of the con- 



tinuum background and misidentified electrons are ex- 
tracted from the data. For all other contributions, Monte 
Carlo (MC) simulation provides the shapes of the dis- 
tributions. The decays B — > D^*'ev are simulated us- 
ing a model based on heavy quark effective theory [l^ . 
The modes B —* D^nev are simulated according to 
the Goity-Roberts model The resonances b — > uev 
heavier than p and u> are implemented according to the 
ISGW2 model 0. Non-resonant b — > uev modes are de- 
scribed by the model of Fazio and Neubert . 

The fit has nine free parameters: B(B° — » p~e + v), 
B(B° — » n~e + v), the normalization of the b — > uev back- 
ground in the two electron-energy ranges (two parame- 
ters), and the normalization of the b — > cev background 
(five parameters, one for each mode). The rates of the 
p , u, and 7T° channel are constrained by the isospin and 
quark model relations T(B° — > p~e + v) = 2T(B + — > 
p°e+v), T(B+ -> p°e + v) = T(B+ -> ue+v), and 
T(B° -> ij-e+v) = 2T(B+ -> ir°e+v). The maximum- 
likelihood fit takes into account the statistical uncertain- 
ties in the on- and off-resonance data and in the proba- 
bility distributions extracted from MC simulations [l5| . 

Projections of the data and fit results for B° — > p~e + v 
are shown in Fig^for the ISGW2 model. A continuum- 
background contribution of 917 ± 73 events in high-i? e 
and 1928 ± 106 in low--E e has been subtracted. Good 
agreement between data and the fit result is seen in each 
of these figures. The fits for the other form- factor cal- 
culations show the same level of agreement. The fit 
quality has been checked with a % 2 test, where bins in 
sparsely populated regions have been combined before 
the x 2 calculation. We obtain \ 2 = 91 f° r 93 degrees 
of freedom for ISGW2, and similarly good fit quality for 
the other form-factor calculations. The signal yields ex- 
tracted from the maximum-likelihood fit in the high-£" e 
region are 321 ± 40 B + — > p°e + v events and 505 ± 63 
B° — » p~e + v events. The resulting branching fractions 
B(B° — > p~e + v) are shown in Fig. EI The five fit param- 
eters describing the b — > cev backgrounds agree with the 
known branching fractions jlpj for B — > Dev, B — > D*ev, 
and B —> D^*'(n)ev within ±9% on average. The two pa- 
rameters describing the size of the background from other 
b — > uev decays agree within 1.5(7 with the predictions 
of the MC simulation. The fit result for the ir modes 
is B(B° -> TT'e+v) = (1.86 ± 0.56 stat .) x 10" 4 for the 
ISGW2 model. 

A summary of all considered systematic uncertainties 
on B(B — > pev) is given in Table llT| The total system- 
atic uncertainty is the quadratic sum of all individual 
ones. Note that the statistical uncertainties in Fig. |21 al- 
ready include the statistical uncertainty in the MC pre- 
dictions. The largest single contribution to the system- 
atic error arises from the uncertainty in the shape of the 
b — * uev background from events other than the signal 
modes. The fraction of b — > uev background events that 
are non-resonant is varied from to 2/3 to estimate this 
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TABLE II: Summary of all contributions to the systematic 
uncertainty on the branching fraction B(B — > pev). 



-2 -1.5 -1 -0.5 0.5 1 1.5 2 
AE [GeV] 



Contribution 


oB p /B p (7o) 


Tracking efficiency 


±5 


Tracking resolution 


±1 


7r efficiency 


±5 


7r energy scale 


±3 


b —* cev background composition 


+1.4, -1.7 


Resonant b —* uev background composition 


+6, —4 


Non-resonant 6 — ► uev background 


±9 


B lifetime 


±1 


Number of BB pairs 


±1.6 


Misidentified electrons 


< ±1 


Electron efficiency 


±2 


B{Y{AS) -> B + B~)/B(T(4S) -> B°B°) 


< ±1 


Isospin and quark model symmetries 


< ±1 


Fit method 


+4,-6 


Total systematic uncertainty 


±14.4 



FIG. 1: Continuum-subtracted data distributions (points 
with error bars) and fit projections (histograms) for M vn (top 
plots) and AE (bottom plots) for the B° — » p~ e + v channel in 
the \ow-E e (left plots) and high-i5 e regions (right plots). The 
fit results are shown for the ISGW2 model. The histograms 
correspond to the true and crossfeed components of the sig- 
nal (open histogram, above and below the dashed line, re- 
spectively) , the background from other b —> uev decays (dark 
shaded region), and b — > cev and other backgrounds (light 
shaded region). 



ISGW2: 

1 1 • 1 1 2.76 + 0.34 + 0.40 

Beyer/Melikhov: 

1 >— ' • 1 1 3.64 ± 0.46 + 0.52 

UKQCD: 

1 1 * 1 1 3.34 + 0.42 ±0.48 

LCSR: 

1 1 • 1 1 3.86 ±0.50 ±0.56 

Ligeti/Wise: 

1 1 • M 2.86 ±0.37 ±0.41 

Combined: 

1 1 1 * 1 1 1 3.29 ± 0.42 ± 0.47 ± 0.60 

2 2.5 3 3.5 4 4.5 

B{B° ^ p-e+v)/W~ 4 

FIG. 2: The B° — » p~e + v branching fraction results us- 
ing five different form-factors calculations. The uncertainties 
shown are statistical, systematic, and (for the combined re- 
sult) theoretical, successively added in quadrature. The com- 
bined result is the average of the five form-factor results. 

uncertainty. The composition of the resonant component 
of other b — ► uev decays has been varied by changing the 
branching fractions for individual resonances by ±50%, 
while keeping the total rate constant. Variations in the 
b — > cev composition contribute much less to the total 
systematic error than variations in the b — > uev compo- 
nent. Possible violations of the isospin and quark model 



constraints are discussed in Refs. |l]| and Their 
contribution to the systematic error is determined by al- 
lowing a ±3% violation. Several fits were performed: 
fitting without the u> mode, without the 7r mode (fixing 
B{B -> 7re^) Q), without the \om-E e region, and with 
different binning. We assign a systematic uncertainty 
for the fit method as half the largest resulting changes 
of the fit result. We have also varied the most impor- 
tant selection requirements and find that the changes in 
B(B — y pev) are consistent with statistical variations as 
determined by a MC simulation. 

A value of \V u b\ is determined by the relation 

\V ub \ = ^B(B° ^ p-e+v)/(f th T B o) , (3) 

where T t h is the predicted form-factor normalization as 
given in Tabled The branching fractions are used sepa- 
rately for each form- factor calculation, as shown in Fig. El 
We use t b o = 1.542 ± 0.016 ps 16] for the B° lifetime. 
The results for \V u b\ are shown in Fig. [3] The combined 
result has been obtained as weighted average of the five 
form- factor results, where the weight is obtained from 
the theoretical uncertainty of each. The theoretical un- 
certainty on the combined result is estimated to be half 
of the full spread of all theoretical uncertainties. 

In conclusion, we have measured the branching fraction 
B(B° -y p~ e+ v) = (3.29 ± 0.42 ± 0.47± 0.60) x 10~ 4 
using isospin constraints and extrapolating to all electron 
energies according to five different form-factor calcula- 
tions. The errors given are statistical, systematic, and 
theoretical, in the order shown. The value of \V u b\ de- 
termined by the same form-factor calculations is \V u b\ = 
(3.64 ± 0.22 ± 0.25 ±° f 6 ) x 10~ 3 . Our results are 
slightly higher (22% for B and 13% for \V ub \) than a pre- 
vious B — y pev result from CLEO [8|, but agree within 
statistical errors. 
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ISGW2: 

3.55 + 0.21 ±0.25 *°™ 
Beyer/Melikhov: 

3.84 ± 0.24 ± 0.27 *°J 
UKQCD: 

3.62 ±0.22 ±0.25 *" 6 
LCSR: 

3.85 ± 0.24 ± 0.27 *° 6 5 7 7 
Lieeti/Wise: 

3.09 ±0.19 ±0.22 *° 4 4 , 2 

Combined: 

3.64 ± 0.22 ± 0.25 



dation. 



FIG. 3: | Viib | determined using five different form-factor cal- 
culations. Only theoretical error bars are shown. The com- 
bined result is also shown at the bottom with statistical, sys- 
tematic, and theoretical uncertainties successively added in 
quadrature. The combined result is the weighted average of 
the five form-factor results, where we have used only the the- 
oretical uncertainties to calculate the weights. 
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